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Abstract 

I discuss the state of the art and outline direction for research in 
event generation for electroweak physics at LEP2 and e + e~ Linear 
Colliders. 

1 Introduction 

With the start of experimentation at LEP1, a new era of precision measure- 
ments in high energy physics had begun. The potential of such machines 
calls for a new level of sophistication in the theoretical tools. The ability to 
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measure parameters of the electroweak standard model at the permille level 
requires precise theoretical predictions of observable quantities. 

Today's sophisticated multi purpose detectors have a complicated geom- 
etry that can not be described by simple acceptance and efficiency functions, 
but require detailed simulation instead. For this purpose, Monte Carlo event 
generators and sophisticated integration algorithms are indispensable. 

The theoretical issues involved in the calculation of differential cross sec- 
tions have been covered by other speakers at this conference |]J and need not 
be repeated here. I will therefore concentrate on the more technical issues 
involved in turning these theoretical predictions into numbers that can be 
compared with experimental results. 

I will start in section Q with reviewing the motivation for Monte Carlo 
methods and I will describe the methods in use today in section ^|, including 
Quasi Monte Carlo in subsection |3.3| . I will continue with radiative correc- 
tions in section f| and discuss specific problems with automated calculations 
of radiative corrections in subsection |4.1| . After reviewing the status of event 
generators for LEP2 in section |^, I describe the new challenges at a Linear 
Collider in section including a discussion of beamstrahlung in section |6~2 . 
After a brief look into the future in section [F] and a few bits and pieces in 
section Bl, I will conclude. 



2 The Need for Monte Carlo 

The increasing energy of e + e~ storage rings and Linear Colliders opens up 
channels with many particles in the final state, that are calculable in the 
electroweak standard model. Continuing the precision tests of the standard 
model at higher energies requires a precise calculation of the cross sections 
for these channels. In addition, these channels provide backgrounds for new 
particle searches, which calls for theoretical control of these cross sections as 
well. 

Numerical integration estimates the integral of a function / by a weighted 
sum of the functions' values sampled on a set {x} of points 

Jdxf{x) = (f)= w ( x )f(*)- (!) 

x£{x} 

Monte Carlo integration is the special case of ([I]), in which {x} is a sam- 
ple of pseudo random numbers and w(x) = l/\{x}\, where \{x}\ denotes 
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the number of points in {x}. This case is of special importance, because 
the law of large numbers guarantees the convergence of (|]) for any /, only 
the scale of the error depends on /. In addition, the statistical nature of 
the method allows reliable error estimates by repeating the evaluation with 
varying random sets and checking the Gaussian distribution of the results. 

Monte Carlo is the only known method that allows the integration of 
differential cross sections for arbitrary final states with arbitrary phase space 
cuts. Monte Carlo event generation is required for realistic simulations of 
the acceptance and efficiency of modern detectors with their complicated 
geometry. Thus, Monte Carlo is the universal tool for turning actions into 
measurable numbers. 

2.1 Discovery vs. Precision Physics 

Discovery physics and precision physics call for different approaches to event 
generation. Precision tests of the standard model require complete calcula- 
tions, that include radiative corrections and "background" diagrams. Since 
the number of diagrams for many particle final states is large (up to 144 for 
four fermions and thousands for six fermions), such calculations are techni- 
cally demanding. Fortunately, the parameter space is restricted and allows 
a comprehensive analysis. 

On the other hand, the search for physics beyond the standard model 
in general does not need radiative corrections other than the initial state 
radiation of photons. This makes the calculations technically easier, but 
the need to cover a lot of models with a vast parameter space creates other 
problems. Since tree level calculations folded with initial state structure 
functions are usually sufficient, computer aided approaches can help to cover 
the models and their parameter spaces. 

3 Twelve Roads from Actions to Answers 

There are twelve different methods for getting answers (cross sections and 
event rates) from actions (the definition of the physical model in perturba- 
tive calculations). These can be factorized in three roads from actions to 
amplitudes and four ways from amplitudes to answers. 
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3.1 The Three Roads from Actions to Amplitudes 



The traditional textbook approach to deriving amplitudes from actions are 
manual calculations, which can be aided by computer algebra tools. The 
time-honored method of calculating the squared amplitudes directly using 
trace techniques is no longer practical for today's multi particle final states 
and has generally been replaced by helicity amplitude methods. Manual 
calculation has the disadvantage of consuming a lot of valuable physicist's 
time, but can provide insights which are hidden by the other approaches 
discussed below. 

The currently most successful and increasingly popular technique is to 
use a well tested library of basic helicity amplitudes for the building blocks 
of Feynman diagrams and to construct the complete amplitude directly in 
the program in the form of function calls. A possible disadvantage is that 
the differential cross section is nowhere available as a formula, but the value 
of such a formula is limited anyway, since they can hardly be printed on a 
single page anymore. 

Automatic calculations are a further step in the same direction. The Feyn- 
man rules (or equivalent prescriptions) are no longer applied manually, but 
encoded algorithmically. This method will become more and more important 
in the future, but more work is needed for the automated construction of ef- 
ficient event generators and for the implementation of radiative corrections 
(see also section Of). 



3.2 The Four Roads from Amplitudes to Answers 

Analytic and semi-analytic calculations have the potential to provide the 
most accurate results. Unfortunately, fully analytic calculations are not fea- 
sible for more than three particles in the final state, that is for most of the 
interesting physics at LEP2 and the Linear Collider. Still, semi-analytic 
calculations, in which some integrations are performed analytically and the 
remaining integrations are done numerically are possible for certain simple 
cuts and provide useful benchmarks with an accuracy unmatched by the 
other methods. 

The most flexible approach is Monte Carlo integration, which converges 
under very general conditions as 
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with N denoting the number of function evaluations. As long as the in- 
tegration variables can be transformed such that the integrand / does not 
fluctuate too wildly, this method can be implemented easily and efficiently 
The quadratic increase on computing resource consumption with the preci- 
sion puts a practical limit on the attainable precision, however. 

Event Generation is a special case of Monte Carlo integration in which 
an ensemble of kinematical configurations is generated that is distributed 
according to the differential cross section. Such ensembles allow the straight- 
forward simulation of the non-perturbative physics of the fragmentation and 
hadronization of strongly interacting particles and of the detector. If the 
weight function is bounded, rejection methods can be used to trivially turn a 
Monte Carlo code into an event generator. Hand tuning is required, however, 
to make the code efficient. 

The rate of convergence in (0) is guaranteed by the law of large numbers 
and can only be improved if we turn away from pseudo random numbers 
and switch to deterministic integration and Quasi Monte Carlo. Empirical 
evidence from known quasi random number sequences suggests that 

(cf. (^D, below) is possible by using point sets that are more uniform than 
both random sets and hyper-cubic lattices in high dimensions. Quasi Monte 
Carlo integration has been applied successfully to four fermion production 
at LEP2 0. Nevertheless, more work is needed, because there are too few 
theorems for realistic function spaces. Phase spaces of varying dimensionality, 
as in branching algorithms, have not been studied at all yet. 



3.3 Quasi Random Numbers 

It is intuitively clear that the best convergence will be gotten by using the 
most uniform point sets. It is less obvious how such uniform point sets look 
like and how to construct them. Let us therefore take a closer look at such 
point sets. More detail and references can be found in ||. 

Let us assume for simplicity that we can map the integration region to 
the n-dimensional unit hypercube: / = [0,1] xn . Using the characteristic 
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function 



x(v\x) = l[e(v*-x*), (4) 

n=i 

we can define the local discrepancy of a point set {x} for each y <E I: 

0d/l{*}) = iAr £ x(yk)-n^- ( 5 ) 

Obviously, the local discrepancy measures how uniformly {x} covers the hy- 
per cube. We can now define various versions of the global discrepancy of the 
point set {x}: 

D m ({x}) = J dy {g{y\{x})) m (6a) 

£) 0O ({ a ;}) = sup|^(y|{ a ;})|. (6b) 
v 

A lower bound for D^, can be established 

DMx}) > C(n) ^ypA (7) 

where C (n) depends only on the dimension of the hypercube and is in partic- 
ular independent of {x}. The global discrepancy D 2 of a regular hypercubic 
lattice and a random point set can easily be calculated 

n 1 

D 2 ({lattice}) = — __ + ... (8a) 

V1 fJ 4-3" I {lattice} | 2 /™ V ; 

D 2 ({random}) = ( l - — 1— - . (8b) 



2 n 3 n J I {random} I ' 

This shows that hypercubic lattices are less uniform than random point sets 
for more than two dimensions. This result is not surprising, however, be- 
cause we know from experience that Monte Carlo works better than uniform 
integration formulae in higher dimensions. 

Unfortunately, while it is intuitively obvious that discrepancy and inte- 
gration error are related, it is much harder to derive mathematically rigorous 
results for realistic integrands, in particular for those with singularities and 
discontinuities. This is the price to pay for potentially faster convergence 
and more research is needed. 
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V k 
Figure 1: Cancellations of real (a) and virtual (b) singularities. 

4 Radiative Corrections 

At high energies, electromagnetic radiative corrections are enhanced by large 
logarithms ln(s/m 2 ) and precision calculations of the hard cross section are 
useless if the radiative corrections are not under control. 

Perturbative calculations at fixed order order in perturbation theory are 
not sufficient, because amplitudes for the emission of photons (figure |l|a) have 
l/tu-soft and l/^-collinear singularities. These singularities are cancelled 
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in the cross section by opposite-sign singularities in the virtual corrections 
(figure |TJb) to the indistinguishable process without emission. While suffi- 
ciently integrated cross sections remain positive, the differential cross section 
without emission must be negative for the cancellation to take effect. Only 
after resummation of the perturbation series, these singularities become in- 
tegrable like I/uj 1-13 and result in strictly positive differential cross sections. 
The combination of today's collider energy and of the energy- and angu- 
lar resolution of today's detectors have made this resummation a practical 
necessity. 

The most popular approach today is the summation of the leading loga- 
rithmic initial state contributions through structure functions 



da 
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an 









D(x + , s)D(x. 



(10) 



The structure functions D can be obtained from finite order perturbation the- 
ory with resummed ("exponentiated") singularities, from Yennie-Frautschi- 
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Suura summation, from the DGLAP evolution equations or from parton 
shower algorithms. 

A drawback of this approach is that the transversal momentum carried 
away by the photons is ambiguous in leading logarithmic approximation. 
For practical purposes, most authors use the l/(p • k) behaviour suggested 
by the most singular pieces. This problem is worrisome for experiments that 
measure the amount of initial state radiation by tagging photons in order to 
select an event sample with effective center-of-mass energy below the nominal 
machine energy. Theoretically most satisfactory are algorithms that match 
the hard matrix element outside the singular region to the resummed cross 
section in the singular region. 

Another universal correction appears in the W + W _ -production close to 
threshold at LEP2. The long range Coulomb interaction of slow W's gives 
rise to a substantial change of the cross section. This correction is easy to 
implement and available in most computer codes. 

The issue of QCD corrections was the subject of some debates during the 
1995 LEP2 workshop. The "naive" QCD correction 



is obtained by integrating over the gluon in the decay of a vector boson into 
quark pairs. The effect on the width and the branching ratios is substantial 
and must be included in some way in the calculations. At the same time, 
this naive approach is theoretically questionable, because the unavoidable 
experimental cuts will invalidate the fully inclusive calculation. Since there 
are again cancellations between real and virtual contributions at work, this is 
not a priori a numerically small effect. Furthermore, the resonant diagrams 
do not form a gauge invariant subset and a more complete analysis is called 
for. 

Such studies are under way || . But until their results are included in the 
standard computer codes, the naive QCD corrections are available in most 
programs because the effect is substantial. 

Weak radiative corrections are certainly more interesting from the physics 
perspective, but they will hardly be measurable at LEP2. Their calculation 
is very tough, because effective Born approximations similar to those suc- 
cessful at LEP1, do not exist for differential cross sections in W + W~ pair 
production. Since complete calculations do not yet exist except for the on- 
shell production of stable W's, gauge dependent ad-hoc resummations of the 




(11) 



8 



propagators are needed. 

For practical purposes at LEP2, the gauge dependence problem has been 
solved by the fermion loop scheme ||, but more complete investigations and 
implementations remain desirable, because stronger effects are expected at 
the Linear Collider. 

4.1 Radiative Corrections in Automatic Calculations 

Radiative corrections in gauge theories pose a particular problem for au- 
tomatic calculations (see [H] for a general review of automatic calculations). 
Formally |7|] , loops can be incorporated in existing tree level generators by re- 
placing the classical action S(A, ip, ijj) with the generating functional T(A, ip, 
of one particle irreducible diagrams. Unfortunately, this approach works only 
for loops consisting of heavy particles. 

As mentioned above, the amplitude for the emission of massless parti- 
cles (figure [IJa) contains infrared and collinear singularities, which are can- 
celled in the cross section by the virtual contributions (figure 0b) for degen- 
erate final states. The problem is that both contributions have to be reg- 
ularized and that the appropriate regularization of the loop in diagram |I]b 
depends on which external leg it is attached to. Therefore it is not possible 
to describe the loop diagrams by a position independent generating func- 
tional F(A,if),'^). More sophisticated algorithms which analyze to topology 
of the graphs are needed. A new approach to this problem for QCD jet cross 
sections has been presented at this conference ||. 

Until this problem is solved, automatic calculations are useful for the 
emission of hard photons which are separated from all charged photons. In- 
tegrated cross sections can be calculated in the leading logarithmic approxi- 
mation only by folding hard cross sections with structure functions (|T0|). 

5 Status of Event Generators for LEP2 

The event generators available for lUIU-physics at LEP2 have been described 
in ||. An updated summary of the properties of the available codes is pre- 
sented in table [I]. 

In ||, the predictions for LEP2 have been compared in detail. In a 
"tuned" comparison with a prescribed calculational scheme, the implemen- 
tations have been tested and the modern dedicated e + e~ — > 4/ codes agreed 
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Table 1: Properties of the available computer codes for H /+ W / ~-physics 
at LEP2. In the 'Type' column 'MC stands for Monte Carlo integration, 
'EG' for event generation, 'SA' for semi-analytic and 'Int.' for determin- 
istic integration. Implemented subsets of diagrams are denoted by 'CC03' 
for doubly resonant W + W~ , 'CC11' for singly resonant W + W~ , 'CC20' for 
final states including electrons or positrons, 'NCnn' for various neutral cur- 
rent diagrams and 'NCC for various neutral and charged current diagrams. 
The implementation of initial state radiation is denoted by 'SF' for struc- 
ture functions, 'FF' for flux functions, 'PS' for parton showers, 'YFS' for 
Yennie-Frautschi-Suura, 'ME' for matrix element and 'O(a)' for one pho- 
ton bremsstrahlung. The 'NQCD' column applies to naive inclusive QCD 
corrections. For fermion masses, '±' denotes massless matrix elements with 
massive kinematics. See || for references and more details. 
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at a level far better than required at LEP2. In a second "unleashed" compar- 
ison, the different theoretical approaches of the codes have been compared 
and the results have shown that the predictions are under control at the level 
required for LEP2. 



6 Towards the Linear Collider 

We have to wait another decade until a Linear Collider will be available for 
experimentation. Nevertheless, the design of the interaction region and of 



detectors needs firm predictions for the expected physics [10] today. 

Regarding two gauge boson physics (which for precision measurements is 
really four fermion physics), most of the LEP2 Monte Carlos in table [3] can 
simply be run at higher energies, provided that they can be interfaced to 
beamstrahlung codes, as discussed in subsection |0] below. 

It is unreasonable to expect deviations of the three gauge boson couplings 
from the Standard Model values (anomalous couplings) that will be measur- 
able at LEP2 (though this assertion has to be checked anyway). This will 
change at the Linear Collider, because reasonable values of (9(l/(167r 2 )) be- 
come accessible and event generators have to support anomalous couplings. 
Fortunately, the majority of the event generators supports them today and 
preliminary cross checks have been satisfactory. 

As mentioned above, weak radiative corrections will also become relevant. 
Here more work is still needed. 



6.1 e + e — ► 6f + n"f 

At a 500 GeV Linear Collider interesting e + e~ — > 6f channels open up. For 
the first time, precision measurements of ti production at a e + e~-collider 
will be possible. An event generator at the signal diagram level, including 
bound state effects, is available. The study of background diagrams will be 
necessary and a cross check of the generator will be desirable. 

At high energies, vector boson scattering (cf. figure |2|), becomes an inter- 
esting e + e~ — > 6f channel as a probe of the electroweak symmetry breaking 
sector. Work in this direction has started this year. 

The general case of e + e~ — > 6/ is a formidable computational challenge, 
because a huge number of diagrams has to be calculated. At the moment 
four approaches are being discussed: 
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Figure 2: VV — > VV scattering. 



1. start from on-shell e + e~ — > VVV, VV — > and e + e~ — > tt Monte 
Carlos and add V —> ff decays in a second step. This approach is 
probably not useful for obtaining a complete calculation in the end, but 
it can provide reasonable descriptions of the most important, resonant 
contributions soon. 

2. extent the EXCALIBUR algorithm for massless four fermion produc- 
tion to six fermions. It is still unclear how to deal with quark masses 
efficiently in this approach. 

3. start with complete calculations of specific final states for which the 
number of diagrams is manageable. Work in this direction has started 
and is showing first promising results [11|. 

4. perform a fully computerized calculation. Work in this direction using 
the GRACE || system has started as well. 

These projects will certainly keep the aficionados of standard model calcula- 
tions entertained for some years. 

6.2 Beamstrahlung 

The experimental environment at the Linear Collider causes a new phe- 
nomenon that event generators have to deal with. The largest (two fermion) 
standard model cross sections are of the order 
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and four and six fermion cross sections are suppressed by further factors of a. 
It is therefore clear that interesting physics at the Linear Collider will require 
large luminosities of the order of 

10 34 cm- 2 sec _1 » lOOflrV 1 , (13) 

where v = 10 7 sec ~ year/7r corresponds to one "effective" year of running. 
These luminosities can only be achieved with extremely dense bunches of 
particles. 

Such bunches experience a strong electromagnetic interaction with non- 
trivial consequences. A desired effect of this interaction for oppositely charged 
particles is the "pinch effect", which increases the luminosity by further colli- 
mating the bunches through the reciprocal attraction. But the same physics 
gives rise to undesirable side effects as well. The deflection of the particles 
in the bunches causes them to loose several percent of their energy as syn- 
chrotron radiation (beamstrahlung). Therefore, the colliding particles will 
have a non-trivial energy spectrum. This spectrum has to be included in 
the event generators for realistic predictions. At the same time, the radi- 
ated photons take part in je ± - and 77-collisions and therefore their energy 
spectrum has to be known as well. 

Quantitatively, the effect of the beamstrahlung is of the same order as 
the effect of initial state radiation. But unlike ordinary bremsstrahlung, 
beamstrahlung can not be treated by ordinary perturbation theory, because 
the underlying physics is the interaction of a particle with all the particles 
in the colliding bunch. Approximate analytical treatments of this collective 
effect exist, but full simulations ||T^| show that they are not adequate. 



The full simulations consume too much computer time and memory for 
directly interfacing them to particle physics Monte Carlos. Also, the input 
parameters collected in table § not familiar to most particle physicists. The 
pragmatical solution of this problem is to describe the result of the simula- 
tions by a simple ansatz, that captures the essential features. 

The approximate solutions motivate a factorized ansatz, which should 
positive and have integrable singularities at the endpoints x e ± — > 1 and x 7 — > 


D PlP2 (xi,x 2 ,s) = d Pl (xi)d P2 (x 2 ) (14) 
d e ±{x) = a S(l -x)+ aia; a2 (l - xf 3 (15) 
d 7 ( x ) = a 4 x a5 (l-x) ae . (16) 
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Table 2: Three prototypical linear collider designs at 500 GeV: SBAND and 
TESLA are DESY's room temperature and superconducting options, XBAND is 
for KEK's and SLAC's plans. 



S-Band, ^fs = 500GeV S-Band, y/s = 500GeV 




1 0.01 lo- 4 io- 6 io- 8 io- 8 io- 6 io- 4 0.01 1 

1 - x e ± x 1 



Figure 3: Quality of the fits for the TESLA design at 500 GeV and 1 TeV. 
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Figure 4: Results of simulations and fits for the three prototypical linear 
collider designs at 500 GeV. 

Figure |3| shows that this ansatz works surprisingly well. It has therefore been 
made available as distribution functions and non-uniform random number 
generators in the circe library |13| . 



The resulting distributions are displayed in figure f|, which highlights the 
substantial differences in the designs. 



7 Monte Carlo Futures 

The craft of Monte Carlo event generation for precision physics at high energy 
e + e~-colliders is practised successfully by many physicists today, as witnessed 
by the many high quality computer codes available for LEP2. Still there are 
areas where technical progress is desirable. 

The most promising direction are generator generators, where amplitudes 
and event generation or integration algorithms are constructed algorithmi- 
cally by a computer program from a Lagrangian for each desired final state. 
The systems available today have come a long way ||, but are still far from 
perfect. 

While the combinatorics of the generation of the amplitude is imple- 
mented in several systems, most of these implementations do not scale well 
to final states with more than four particles. The number of diagrams will 
scale like n\ with the number n of particles. For small n, the resource con- 



15 



sumption of this factorial growth can be matched by faster hardware in the 
future, but this is not obvious for the numerical stability of the code. Gauge 
theories have a good high energy behaviour because of strong cancellations 
among individual diagrams. The numerical problems resulting from these 
cancellations have to be controlled. This is only possible with more sophis- 
ticated algorithms that regroup and partially factorize the amplitudes. Such 
algorithms are non trivial and will need further research. Here progress is 
particularly important for radiative corrections. 

A second area of research are adaptive methods for event generation and 
integration. Except for special cases, where the singularities in the ampli- 
tudes are known in advance, today's systems still need human intervention 
for finding optimal phase space variables that minimize the sampling error. 
Here more sophisticated algorithms are needed as well. 

Technical advances in this direction will allow to efficiently produce reli- 
able Monte Carlo codes for precision physics at high energy Linear Colliders. 
Hopefully, this will also propel the state of the art of event generation for 
discovery physics to the same level. 

Another promising direction for research is the investigation of Quasi 
Monte Carlo methods for event generation and integration. Here more expe- 
rience with realistic applications is needed. 

There are however some areas where progress will be slow. In particular 
the interface of the perturbative precision calculations with the non pertur- 
bative simulation of strongly interacting final states is poorly understood. 
A lot of the sophistication in the calculation of interfering contributions is 
lost when the perturbative amplitude is matched to the classical simulation 
of the fragmentation and hadronization process. At LEP2, the problems of 
color reconnection and Bose-Einstein correlations are the limiting factor for 
the H^-mass measurements. When LEP2 data become plentiful, they might 
help to improve phenomenological models and give a better control of the 
systematical error. 

8 Bits and Pieces 

Before concluding, I want to take this opportunity to advertize a welcome 
addition to the family of pseudo random number generators available for 
Monte Carlo integration and event generation. 

It is well known that bad random number generators will spoil any simula- 
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tion, while (some) good random number generators can consume macroscopic 
fractions of the total computer time. 



Recently, Donald Knuth made errata [13] for his textbook available, which 
contain a gem of a portable generator. It is an extremely fast implementation 
of a lagged Fibonacci generator Xj = (Xj^ wo — X,_ 37 ) mod 2 30 , which 
is portable even for systems which, like FORTRAN, offer no unsigned 32-bit 
arithmetic. The generator passes all statistical tests, even (with a slight 
speed penalty) the birthday spacings test. But the most interesting property 
is the innovative initialization algorithm for which one can prove that it 
will generate 2 30 — 2 statistically independent sequences from a simple 30- 
bit integer seed. For most other generators, the statistical independence for 
different seeds is much harder to control. 

The practical consequence is that Knuth's new algorithm allows parallel 
execution and reliable Monte Carlo estimation of errors, without having to 
worry about the statistical independence of the generated samples. 

9 Conclusions 

The Monte Carlo codes for precision physics at LEP2 physics are in excellent 
shape and the craft of their construction is well understood. Tools for the 
design of the Linear Collider are already available, but there is work left to 
do to make them as comprehensive as the tools available for LEP2. 

Technical progress is still needed for streamlining the computational ap- 
proaches that will eventually allow us to investigate the nature of electroweak 
symmetry breaking at the Linear Collider. 
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